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Abstract 

Video content delivery over wireless networks is expected to grow drastically in the coming years. 
In this paper, we investigate the challenging problem of video over cognitive radio (CR) networks. 
Although having high potential, this problem brings about a new level of technical challenges. After 
reviewing related work, we first address the problem of video over infrastructure-based CR networks, 
' and then extend the problem to video over non-infrastructure-based ad hoc CR networks. We present 

fSJ formulations of cross-layer optimization problems as well as effective algorithms to solving the problems. 

m 

The proposed algorithms are analyzed with respect to their optimality and validate with simulations. 
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I. Introduction 

Video content delivery over wireless networks is expected to grow drastically in the coming 



years. The compelling need for ubiquitous video content access will significantly stress the 
capacity of existing and future wireless networks. To meet this critical demand, the Cognitive 
Radio (CR) technology provides an effective solution that can effectively exploit co-deployed 
networks and aggregate underutilized spectrum for future video-aware wireless networks. 

The high potential of CRs has attracted substantial interest. The mainstream CR research 
has focused on developing effective spectrum sensing and access techniques (eg., see 0], [|2l). 
Although considerable advances have been achieved, the important problem of guaranteeing 
application performance has not been well studied. We find video streaming can make excellent 
use of the enhanced spectrum efficiency in CR networks. Unlike data, where each bit should be 
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delivered, video is loss-tolerant and rate-adaptive O, @|. They are highly suited for CR net- 
works, where the available bandwidth depends on primary user transmission behavior. Graceful 
degradation of video quality can be achieved as spectrum opportunities evolve over time. 

CR is an evolving concept with various network models and levels of cognitive functional- 
ity JD. ED- IEEE 802.22 Wireless Regional Area Networks (WRAN) is the first CR standard for 
reforming broadcast TV bands, where a base station (BS) controls medium access for customer- 
premises equipments (CPEs) [5]. Therefore, we first consider multicasting scalable videos in such 
an infrastructure-based CR network. The spectrum consists of multiple channels, each allocated 
to a primary network. The CR network is co-located with the primary networks, where a CR 
BS seeks spectrum opportunities for multicasting multiple video streams, each to a group of 
secondary subscribers. The problem is to exploit spectrum opportunities for minimizing video 
distortion, while keeping the collision rate with primary users below a prescribed threshold. We 
consider scalable video coding, such as fine-grained-scalability (FGS) and medium grain scalable 
(MGS) videos J6]|, 0. We model the problem of CR video multicast over the licensed channels as 
a mixed integer nonlinear programming (MINLP) problem, and then develop a sequential fixing 
algorithm and a greedy algorithm to solve the MINLP, while the latter has a low computational 
complexity and a proved optimality gap j8]|. 

We then tackle the problem of video over multi-hop CR networks, e.g., a wireless mesh 
network with CR-enabled nodes. This problem is more challenging than the problem above 
due to the lack of infrastructure support. We assume each secondary user is equipped with two 
transceivers. To model and guarantee end-to-end video performance, we adopt the amplify-and- 
forward approach for video data transmission, which is well-studied in the context of cooperative 
communications [9]. This is equivalent to setting up a "virtual tunnel" through a multi-hop multi- 
channel path. The challenging problem, however, is how to set up the virtual tunnels, while 
the available channels at each relay evolve over time due to primary user transmissions. The 
formulated MINLP problem is first solved using a centralized sequential fixing algorithm, which 
provides upper and lower bounds for the achievable video quality. We then apply dual decom- 
position to develop a distributed algorithm and prove its optimality as well as the convergence 
condition ifTOll . 

The rest of the paper is organized as follows. We review related work in Section [TT] and 
present preliminaries in Section Hill We examine video over infrastructure-based CR networks in 
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Section [IV] and over multi-hop CR networks in Section [V] We concludes the paper in Section [VI] 
with a discussion of open problems. 

II. Background and Related Work 

The high potential of CRs has attracted considerable interest form both industry, government 
and academia HI, ifTTI . The mainstream CR research has been focused on spectrum sensing 
and dynamic spectrum access issues. For example, the impact of spectrum sensing errors on the 
design of spectrum access schemes has been addressed in several papers |[T2 l -|fT4 | . The approach 
of iteratively sensing a selected subset of available channels has been developed in the design 
of CR MAC protocols [fT31 . |[T6ll . The optimal trade-off between the two kinds of sensing errors 
is investigated comprehensively and addressed in depth in 021. 

The important issue of QoS provisioning in CR networks has been studied only in a few 
papers lfl"5l . [fTTl . where the objective is still focused on the so-called "network-centric" metrics 
such as maximum throughput and delay lfl"3~1 . lfl"5l . In lfl"3~1 . an interesting delay throughput trade- 
off for a multi-cell cognitive radio network is derived, while the goal of primary user protection 
is achieved by stabilizing a virtual "collision queue". In ifTTll . a game-theoretic framework 
is described for resource allocation for multimedia transmissions in spectrum agile wireless 
networks. In this interesting work, each wireless station participates in a resource management 
game, which is coordinated by a network moderator. A mechanism-based resource management 
scheme determines the amount of transmission opportunities to be allocated to various users on 
different frequency bands such that certain global system metrics are optimized. 

The problem of video over CR networks has been addressed only in a few recent papers. 
In |[T8l . a priority virtual queue model is adopted for wireless CR users to select channel and 
maximize video qualities. In |[T9ll , the impact of system parameters residing in different network 
layers are jointly considered to achieve the best possible video quality for CR users. The problem 
is formulated as a Mini-Max problem and solved with a dynamic programming approach. In ll20ll . 
Ali and Yu jointly optimize video parameter with spectrum sensing and access strategy. A rate- 
distortion model is adopted to optimize the intra-mode selection and source-channel rate with 
a partially observable Markov decision process (POMDP) formulation. In [21J, video encoding 
rate, power control, relay selection and channel allocation are jointly considered for video over 
cooperative CR networks. The problem is formulated as a mixed-integer nonlinear problem and 
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solved by a solution algorithm based on a combination of the branch and bound framework and 
convex relaxation techniques. 

Video multicast, as one of the most important multimedia services, has attracted considerable 
interest from the research community. Layered video multicast has been researched in the 
mobile ad hoc networks ll22ll . Il23l and infrastructure-based wireless networks 10, [|24l . A greedy 
algorithm is presented in |[24]| for layered video multicast in WiMAX networks with a proven 
optimality gap. 

A few recent works Il25l - ll27l have studied multi-hop CR networks. The authors formulate 
cross-layer optimization problem considering factors from the PHY up to the transport layer. The 
dual decomposition technique [|28l . ||29l is adopted to develop distributed algorithm. We choose 
similar methodology in our work and apply it to the more challenging problem of real-time 
video streaming. 

III. System Model and Preliminaries 

A. Primary Network 

We consider a spectrum band consisting of M orthogonal channels with identical band- 
width ll30ll . OTl . We assume that the M channels are allocated to K primary networks, which 
cover different service areas. A primary network can use any of the M channels without 
interfering with other primary networks. We further assume that the primary systems use a 
synchronous slot structure as in prior work JTJ, |fT5l . Due to primary user transmissions, the 
occupancy of each channel evolves following a discrete-time Markov process, as validated by 
recent measurement studies HI, lfl"5l . Il32ll . 

In primary network k, the status of channel m in time slot t is denoted by S^ n {t) with idle 
(i.e., S^Jt) = 0) and busy (i.e., S^(t) = 1) states. Let X^ L and n k m be the transition probability 
of remaining in state and that from state 1 to 0, respectively, for channel m in primary network 
k. The utilization of channel m in primary network k, denoted by 77^ = Pr(S^ = 1), is 

Note that in infrastructure-based CR networks and cooperative CR networks, we assume there 
is only one K — 1 primary network. In infrastructure-based CR networks introduced in the 
section [TV] we adopt iV as the number of licensed channels since M is denoted as the number 
of modulation-coding schemes. 



Spec 
i 

..- 

1 

1 

i 
\ 

\ 

\ 

y 

/ 
\ 

Spec 


:trum Spectrum 
Primary Network 1 (Channel 1) ^ ) 

^ \" 

x f\ Primary Network N (Channel N) ) 

<$f „„.^--' 

r - .— — w 
;trum Spectrum 




^Primary M "W SSL,* ^^CR 
^network A P^ ork M P^ ork ft) network 
s b f a S f on ^ s b t a ^ on H user 



Fig. 1. An infrastructure-based CR network collocated with N primary networks. 

B. Infrastructure-based CR Networks 

As shown in Fig. \T\ we consider a CR base station multicasts G real-time videos to G 
multicast groups, each of which have N g users, g — 1, 2, • • • , G. The base station seeks spectrum 
opportunities in the N channels to serve CR users. In each time slot t, the base station selects a 
set of channels Ai(t) to sense and a set of channels Aiif) to access. Without loss of generality, 
the base station has |^4i(t)| transceivers such that it can sense channels simultaneously. 

Note that a time slot and channel combination, termed a tile, is the minimum unit for resource 
allocation. 

We adopt the same time-slot structure as in JT), lfl6ll - , which is illustrated in Fig. [2l At the 
beginning of each time slot, the base station senses channels in Ai(t) and then chooses a set of 
available channels for opportunistic transmissions based on sensing results. After a successful 
transmission, the base station will receive an ACK from the user with the highest SNR in the 
target multicast group. Without loss of generality, we assume that each CR network user can 
access all the available channels with the channel bonding/aggregation techniques [|33l . Il34ll . 
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Fig. 2. The structure of a time slot. 




Fig. 3. Illustration of the multi-hop video CR network architecture. 

C. Multi-hop CR networks 

As shown in Fig. [3l we also consider a multi-hop CR network that is co-located with the 
primary networks, within which S real-time videos are streamed among N CR nodes. Let U k 
denote the set of CR nodes that are located within the coverage of primary network k. A video 
session / may be relayed by multiple CR nodes if source z% is not a one-hop neighbor of 
destination d[. We assume a common control channel for the CR network |fT511 . We also assume 
the timescale of the primary channel process (or, the time slot durations) is much larger than 
the broadcast delays on the control channel, such that feedbacks of channel information can be 
received at the source nodes in a timely manner. 

The time slot structure is the same as that in infrastructure-based CR networks. In the sensing 
phase, one transceiver of a CR node is used to sense one of the M channels, while the other is 
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Fig. 4. The cut-through switching model for video data. 



tuned to the control channel to exchange channel information with other CR users. Each video 
source computes the optimal path selection and channel scheduling based on sensing results. In 
the transmission phase, the channels assigned to a video session / at each link along the path 
form a virtual "tunnel" connecting source z\ and destination d\. As illustrated in Fig. 0] each 
node can use one or more than one channels to communicate with other nodes using the channel 
bonding/aggregation techniques ll33l . P4ll . When multiple channels are available on all the links 
along a path, multiple tunnels can be established and used simultaneously for a video session. 
In the acknowledgment phase, the destination sends ACK to the source for successfully received 
video packets through the same tunnel. 

We adopt amplify-and-forward for video transmission Q- During the transmission phase, one 
transceiver of the relay node receives video data from the upstream node on one channel, while 
the other transceiver of the relay node amplifies and forwards the data to the downstream node 
on a different, orthogonal channel. There is no need to store video packets at the relay nodes. 
Error detection/correction will be performed at the destination node. As a result, we can transmit 
through the tunnel a block of video data with minimum delay and jitter in one time slot. 

D. Spectrum Sensing 

Two types of sensing errors may occur during the sensing process. A false alarm may lead to 
waste a spectrum opportunity and a miss detection may causes collision with primary users. In 
a multi-hop CR network, the sensing results from various users may be different. Denote H as 



7 



the hypothesis that channel to in primary network k is idle, and Hi the hypothesis that channel 
to in primary network k is busy in time slot t. The conditional probability that channel to is 
available in primary network k, denoted by a^(t), can be derived as, 



a*(t) = Pr(H \Wt = er, ielCM 

\u k m \-u>L Pr(^iK 



(0*)' 



(2) 



Pr(# |7T*, 

where 6>™ represents a specific sensing result (0 or 1), is the subset of users in U k (i.e., the 
set of CR nodes that are located within the coverage of primary network k) that sense channel 
m, u k m is the number of users in observing channel to is idle, 7r^ represents the history of 
channel to in primary network k, and ip k m and 0^ are defined as: 

<A = pZ^ThI = when AT = 

0™ - P(wr=llHo) - wnen y. - i. 

Based on the Markov chain channel model, we have ((U), which can be recursively expanded: 

Pr(floK) = A£,c&(t - 1) + /4 [l - a ^(i - 1)] 
Pr(ff 1 |7r*) = l-Pr(ff |O. 

£. Video Performance Measure 

Both FGS and MGS videos are highly suited for dynamic CR networks. With FGS or MGS 
coding, each video I is encoded into one base layer with rate R\ and one enhancement layer 
with rate Rf . The total bit rate for video I is Ri = R\ + Rf. 

We consider peak-signal-noise-ratio (PSNR) (in dB) of reconstructed videos. As in prior 
work [6], the average PSNR of video /, denoted as Q u can be estimated as: 

Qi(Ri) = Q\ + fo{Ri - R\) = Q° l + 0iRi, (5) 

where Q\ is the resulting PSNR when the base layer is decoded alone, $ a constant depending 
on the video sequence and codec setting, and Q° = Q\ — f3iR\. We verified the model © with 
several test video sequences using the MPEG-4 FGS codec and the H.264/SVC MGS codec and 
found it is highly accurate. 

Due to the real-time nature, we assume that each group of pictures (GOP) must be delivered 
during the next GOP window, which consists of N G time slots. Beyond that, overdue data from 
the current GOP will be useless and will be discarded. In infrastructure-based network, G video 

8 



stream are multicast to G groups of CR user, so we choose the group index g instead of video 
session index /. 



In this section, we examine the problem of video over infrastructure-based CR networks. We 
consider cross-layer design factors such as scalable video coding, spectrum sensing, opportunistic 
spectrum access, primary user protection, scheduling, error control and modulation. We propose 
efficient optimization and scheduling algorithms for highly competitive solutions, and prove the 
complexity and optimality bound of the proposed greedy algorithm. 

A. Network Model 

1) Spectrum Access : At the beginning of each time slot t, the CR BS senses the M 
channels and compute a n (t) for each channel n. Based on spectrum sensing results, the base 
station determines which channels to access for video streaming. We adopt an opportunistic 
spectrum access approach, aiming to exploit unused spectrum while probabilistically bounding 
the interference to primary users. 

Let 7„ G (0, 1) be the maximum allowed collision probability with primary users on channel n, 
and Pn(t) the transmission probability on channel n for the base station in time slot t. The proba- 
bility of collision caused by the base station should be kept below 7„, i.e., p^(t) [1 — a n (t)] < 7„. 
In addition to primary user protection, another important objective is to exploit unused spectrum 
as much as possible. The transmission probability can be determined by jointly considering both 
objectives, as 



If Pn(t) = 1, channel n will be accessed deterministically. If Pn(t) = 7„/[l — a n (t)] < 1, channel 
n will be accessed opportunistically with probability Pn(t). 

2) Modulation-Coding Schemes : At the PHY layer, we consider various modulation and 
channel coding combination schemes. Without loss of generality, we assume several choices of 
modulation schemes, such as QPSK, 16-QAM and 64-QAM, combined with several choices of 
forward error correction (FEC) schemes, e.g., with rates 1/2, 2/3, and 3/4. We consider M unique 



IV. Video over Infrastructure Based CR Networks 




(6) 
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combinations of modulation and FEC schemes, termed Modulation-Coding (MC) schemes, in 
this paper. 

Under the same channel condition, different MC schemes will achieve different data rates and 
symbol error rates. Adaptive modulation and channel coding allow us to exploit user channel 
variations to maximize video data rate under a given residual bit error rate constraint. When a 
user has a good channel, it should adopt an MC scheme that can support a higher data rate. 
Conversely, it should adopt a low-rate MC scheme when the channel condition is poor. Let 
{MC m } m= i r .M be the list of available MC schemes indexed according to their data rates in the 
increasing order. We assume slow fading channels with coherence time larger than a time slot. 
Each CR user measures its own channel and feedbacks measurements to the base station when 
its channel quality changes. At the beginning of a time slot, the base station is able to collect 
the number n giTn of users in each multicast group g who can successfully decode MC m signals 
for m = 1, 2, ■ ■ ■ , M. 

Since the base layer carries the most important data, the most reliable MC scheme MCb( g ) 
should be used, where b(g) = maxj{z : n g> i = N g }, for all g. Without loss of generality, we 
assume that the base layer is always transmitted using MC\. If a user's channel is so poor that 
it cannot decode the MC\ signal, we consider it disconnected from the CR network. We further 
divide the enhancement layer into M sub-layers, where sub-layer m has rate R e gtm and uses 
MC m . Assuming that MC m can carry 6 9 m bits of video g in one tile, we denote the number 
of tiles for sub-layer m of video g as l g , m > 0. We have 

M M 

R g — ^ ] Rg,m = ^ ] bg,mlg,m- (7) 
m=l m=l 

3) Proportional Fair Allocation : Since we consider video quality in this paper, we define 
the utility for user i in group g as U 9yi = \ogQ 3)i = log (Q b g + /3 g R e g (i)), where R e g {i) is the 
received enhancement layer rate of user i in group g. 

The total utility for group g is U g — X]i=i ^g,i- Intuitively, a lower layer should use a lower 
(i.e., more reliable) MC scheme. This is because if a lower layer is lost, a higher layer cannot 
be used at the decoder even if it is correctly received. Considering the user classification based 
on their MC schemes, we can rewrite U g as follows ll24l : 

M / k \ 

Ug = - n **+0 l0 § Wo + & E Km , (8) 

k=l \ m=l J 
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where n g M+i = 0. The utility function of the entire CR video multicast system is 

G 

U = Y,U 9 . (9) 

9=1 

Maximizing U will achieve proportional fairness among the video sessions ll35l 

B. Optimized Video Multicast in CR Networks 

1 ) Outline of the Proposed Approach : As discussed, the CR video multicast problem is 
highly challenging since a lot of design choices are tightly coupled. First, as users see different 
channels, such heterogeneity should be accommodated so that a user can receive a video quality 
commensurate to its channel quality. Second, we need to determine the video rates before 
transmission, which, however, depend on future channel evolution and choice of MC schemes. 
Third, the trade-off between primary user protection and spectrum utilization should guide the 
scheduling of video packets to channels. Finally, all the optimization decisions should be made 
in real-time. Low-complexity, but efficient algorithms are needed, while theoretical optimality 
bounds would be highly appealing. 

To address heterogeneous user channels, we adopt FGS to produce a base layer with rate R b g 
and an enhancement layer with rate R g . Without loss of generality, we assume R g is prescribed for 
an acceptable video quality, while R g is set to a large value that is allowed by the codec. During 
transmission, we determine the effective rate for each enhancement layer R e < R g depending 
on channel availability, sensing, and MC schemes Q The optimal partition of the enhancement 
layer should be determined such that each sub-layer uses a different MC scheme. 

We determine the optimal partition of enhancement layers, the choices of MC schemes, and 
video packet scheduling as follows. First, we solve the optimal partition problem for every GoP 
based on an estimated (i.e., average) number of available tiles T e in the next GoP window that 
can be used for the enhancement layer, using algorithm GRD1 with complexity 0(MGT e ). 
The tile allocations are then dynamically adjusted in each time slot according to more recent 
(and thus more accurate) channel status using algorithm GRD2, with complexity O(MGK), 
where K <C T e . Second, during each time slot, video packets are scheduled to the available 
channels such that the overall system utility is maximized. The TSA algorithm has complexity 

'The proposed approach can also be used for streaming stored FGS video. 
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0(N\ogN). Both GRD2 and TSA have low complexity and are suitable for execution in each 
time slot. 

In real-time video, overdue packets generally do not contribute to improving the received 
quality. We assume that the data from a GoP should be be delivered in the next GoP window 
consisting of T GoP time slotso Since the base layer is essential for decoding a video, we assume 
that the base layers of all the videos are coded using MC\. For the M sub-layers of the 
enhancement layer, a more important sub-layer will be coded using a more reliable (i.e., lower 
rate) MC scheme. At the beginning of each GoP window, all the base layers are transmitted 
using the available tiles. Retransmissions will be scheduled if no ACK is received for a base 
layer packet. After the base layers are transmitted, we allocate the remaining available tiles in 
the GoP window for the enhancement layer. The same rule applies to the enhancement sub- 
layers, such that a higher sub-layer will be transmitted if and only if all the lower sub-layers 
are acknowledged. This is due to the decoding dependency of layered video. 

In each time slot t, the base station opportunistically access every channel n with probability 
Pn(t) given in ©. Specifically, for each channel n, the base station generates a random number 
x n (t), which is independent of the channel history 9 n (t) and uniformly distributed in [0,1]. If 
%n(t) < Pn{t)> me most important packet among those not ACKed in the previous GoP will 
be transmitted on channel n. If an ACK is received for this packet at the end of time slot t, 
this packet is successfully received by at least one of the users and will be removed from the 
transmission buffer. Otherwise, there is a collision with primary user and this packet will remain 
in the transmission buffer and will be retransmitted. 

In the following, we describe in detail the three algorithms. 

2) Enhancement Layer Partitioning and Tile Allocation : As a first step, we need to determine 
the effective rate for each enhancement layer R e < R e g . We also need to determine the optimal 
partition of each enhancement layer. Clearly, the solutions will be highly dependent on the 
channel availability processes and sensing results. 

Recall that the base layers are transmitted using MC\ first in each GoP window. The remaining 
available tiles can then be allocated to the enhancement layers. We assume that the number of 
tiles used for the enhancement layers in a GoP window, T e , is known at the beginning of the 

2 The proposed approach also works for the more general delay requirements that are multiple GoP windows. 
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GoP window. For example, we can estimate T e by computing the total average "idle" intervals 
of all the N channels based on the channel model, decreased by the number of tiles used for 
the base layers (i.e., R b g /b g ^). We then split the enhancement layer of each video g into M 
sub-layers, each occupying l g ^ m tiles when coded with MC m , m — 1, 2, • • • , M. 

Letting / = [/ x 1; ■ • • , h,M> ^2,1, • • • Ig,m] denote the tile allocation vector, we formulate an 
optimization problem OPT-Part as follows. 

G M 

maximize: U{1) = ^2^2(n g:k - n fljfc+ i) x log 

3=1 fc=l 

G M 

subject to: ^ ^ h,™ < T e (11) 

3=1 m=l 

M 

^ < /.';. ge[l,---,M] (12) 

m=l 

l g ,m>0, m e [l,---,M],ge [1,---,G]. (13) 

OPT-Part is solved at the beginning of each GoP window to determine the optimal partition 
of the enhancement layer. The objective is to maximize the overall system utility by choosing 
optimal values for the i s>m 's. We can derive the effective video rates as R e = Y^m=i bg,m,lg,m- The 
formulated problem is a MINLP problem, which is NP-hard [|24l|. In the following, we present 
two algorithms for computing near-optimal solutions to problem OPT-Part: (i) a sequential fixing 
(SF) algorithm based on a linear relaxation of (flOl) . and (ii) a greedy algorithm GRD1 with proven 
optimality gap. 

a) A Sequential Fixing Algorithm : With this algorithm, the original MINLP is first lin- 
earized to obtain a linear programming (LP) relaxation. Then we iteratively solve the LP, while 
fixing one integer variable in every iteration ll26ll . Il36l . We use the Reformulation-Linearization 
Technique (RLT) to obtain the LP relaxation [37 J. RLT is a technique that can be used to produce 
LP relaxations for a nonlinear, nonconvex polynomial programming problem. This relaxation 
will provide a tight upper bound for a maximization problem. Specifically, we linearize the 
logarithm function in (fTOl) over some suitable, tightly-bounded interval using a polyhedral outer 
approximation comprised of a convex envelope in concert with several tangential supports. We 
further relax the integer constraints, i.e., allowing the Z fl>m 's to take fractional values. Then we 
obtain an upper-bounding LP relaxation that can be solved in polynomial time. Due to lack of 



Qg + fig bg,mlg,r 



m=l 



(10) 
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TABLE I 

The Sequential Fixing (SF) Algorithm 



Use RLT to linearize the original problem 



2 



Solved the LP relaxation 



3 



Suppose Ig.m is the integer variable with the minimum 
(Rs.aI ^ lg,m) or {l g ^ - \ l g ,m\) value among all Z 9 , m 
variables that remain to be fixed, round it up or down to the 



nearest integer 



4: 
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If all Ig.m's are fixed, got to Step 6 

Otherwise, reformulate a new relaxed LP with the newly 



fixed / 



g,m 



variables, and go to Step 2 



6: 



Output all fixed l g , m variables and R e g — Ylm—i bg,mlg, 



space, we refer interested readers to ll37l for a detailed description of the technique. 

We next solve the LP relaxation iteratively. During each iteration, we find the l g ^ which 
has the minimum value for (|7§,ml — lg,m) or (/§,m — Us,mJ) among all fractional / g , m 's, and 
round it up or down to the nearest integer. We next reformulate and solve a new LP with ^„ 
fixed. This procedure repeats until all the / s m 's are fixed. The complete SF algorithm is given in 
Table HI The complexity of SF depends on the specific LP algorithm (e.g., the simplex method 
with polynomial-time average-case complexity). 

b) A Greedy Algorithm : Although SF can compute a near-optimal solution in polynomial 
time, it does not provide any guarantee on the optimality of the solution. In the following, we 
describe a greedy algorithm, termed GRD1, which exploits the inherent priority structure of 
layered video and MC schemes and has a proven optimality bound. 

The complete greedy algorithm is given in Table HH where R = J2 g=1 R e g is the total rate of 
all the enhancement layers and e*j is a unit vector with "1" at the i-th location and "0" at all 
other locations. In GRD1, all the Z 5im 's are initially set to 0. During each iteration, one tile is 
allocated to the m-th sub-layer of video g. In Step 4, lm,g is chosen to be the one that achieves 
the largest increase in terms of the "normalized" utility (i.e., [U(l + e gtm ) — U(i)]/[b gim + R/T e )) 
if it is assigned with an additional tile. Lines 6, 7, and 8 check if the assigned rate exceeds the 
maximum rate R e . GRD1 terminates when either all the available tiles are used or when all the 
video data are allocated with tiles. In the latter case, all the videos are transmitted at full rates. 
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TABLE II 
The Greedy Algorithm (GRD1) 



1: 

2: Initialize A = {1,2, ••• ,G} 

3: WHILE (j2% Em=i k,™ < T e and A is not empty) 

4: Find l gt7 n that can be increased by one: 

/ U(T+e g , m )~U(l) \ 

ei,,rh = argmax 96 A,m€[i,-'-,A/] | — bg ^ +R /Te — J 

5. I — I ~\~ Cg,rh 

6: IF (Em h,™h, m > />•:,) 

7. / — / &g,m 

8: Delete g from A 

9: END IF 

10: END WHILE 



We have the following Theorem for GRD1. 

Theorem 1: The greedy algorithm GRD1 shown in Table [n] has a complexity 0(MGT e ). It 
guarantees a solution that is within a factor of (1 — e -1 / 2 ) of the global optimal solution. 

Proof: (i) Complexity: In Step 4 in Table UH it takes O(MG) to solve for e g< rh- Since each 
iteration assigns one tile to sub-layer m of group g, it takes T e iterations to allocate all the 
available tiles in a GoP window. Therefore, the overall complexity of GRD1 is 0(MGT e ). 

(ii) Optimality Bound: This proof is extended from a result first shown in [[241 for layered 
videos. We first show a property of group utility U g (l), which will be used in the proof of the 
optimality gap. For two vectors and P g , letting A = U g (lg) — U g {fi g ), we have 



A 



M 



n g>k+ i) x log 1 + 



fc=i 

M k 



Y2m=l 0gbg,m(l gt m ^g,m) 
Qg + Em=l Pgbg,mlg t , 



g,m 



k=l m=l 

M M 



2 ) 

g,m/ 



g,mJ 



{n g , k ~ n gM+ i) x log \1 + 8 g b g>m / 
{n g , k - %,fe+i) x log 1 + /3 g b g , m / 



k=l m=l 
M 

m=l 



n b + V 8 b I 2 

V 9 / 4 l J g u g,m L g ^. 

m=l 
k 

Q g + ^ ] figbg,mlg t r. 



m=l 



(14) 



where y + = max{0,?/}. The first inequality is due to the concavity of logarithm functions. 
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Next we prove the optimality bound. Let l t be the output of GRD1 after t iterations. Let the 
utility gap between the optimal solution and the GRD1 solution be F t = U(l*) — U(l t ), and 
£g,m{t) the argument found in Step 4 of GRD1 after t iterations. We have l t = l t _i + e^(i) and 

< E £(£m - W^S-i + e g , m {t)) - [/(d)] 

g m 

g m 9> m \ I / e 

[/(ft - utft-i) 



< 



The first inequality is due to (fl4l) and the second inequality follows Step 4 of GRD1. It follows 

m that J2 g E m < T e and £ 9 E m «,« < R - We have < - ^ w"^. - 
Solving for F t , we have F t < F t ^ {1 - [& 5iA (t) + R/T e ] /(2R)}. 

Suppose the WHILE loop in Table |n] has been executed k times when the solution is obtained. 

F k < F k ^{l-[b^(k) + R/T e ]/(2R)} 

k 

< F J[{l-[b Sifil (t)+R/T e ]/(2R)} 

t=i 

{k 
l-l/(2kR)J2[b Stlh (t) + R/T e 
t=i 

The WHILE loop exits when one or both of two constraints are violated. If ^ J2 m / g m < T e is 
violated, there is no tile that can be used. Therefore k > T e and Ylt=i R/T e > R- If constraint 
11 A is not empty" is violated, all the videos have been allocated sufficient number of tiles and 
will be transmitted at full rates. We have Ylt=i bg,rh(t) > R m this case. It follows that 

7, 

{k 
l-l/(2kR)J2[b^(t)+R/T e 
t=i 

< F [1 - 1/(2*)]* < Foe' 1 / 2 . 

Since F = U(l*), we have U(T k ) > (1 - e~ 1/2 )?7(f*). Therefore, we conclude that the GRD1 
solution is bounded by (1 - e^^Uil*) and U{1*). ■ 



16 



c) A Refined Greedy Algorithm : GRD1 computes lg jTn 's based on an estimate of network 
status S(t) in the next T GoP time slots. Due to channel dynamics, the computed Z 9)m 's may not 
be exactly accurate, especially when T GoP is large. We next present a refined greedy algorithm, 
termed GRD2, which adjusts the Z gi?Ti 's based on more accurate estimation of the channel status. 

GRD2 is executed at the beginning of every time slot. It estimates the number of available 
tiles T e {t) in the next T est successive time slots, where 1 < T est < T Go p is a design parameter 
depending on the coherence time of the channels. Such estimates are more accurate than that in 
GRD1 since they are based on recently received ACKs and recent sensing results. Specifically, 
we estimate T e (t) using the belief vector a{t) in time slot t. Recall that a n (t)'s are computed 
based on the channel model, feedback, sensing results, and sensing errors, as given in (0, and 
©. For the next time slot, a n (t + 1) can be estimated as a n (t + 1) = X n a n (t) + fi n [l — a n {t)} = 
(A n — /i n )a„(i) + fi n . Recursively, we can derive a n (t + r) for the next r time slots. 

a n (t + r) = (A n - i2 n ) T a n (t) + fi n -. — y^T — ^T- ( 15 ) 

1 - (A n - /I n ) 

At the beginning section of a GoP window, all the base layers will be firstly transmitted. 
We start the estimation after all the base layers have been successfully received (possibly with 
retransmissions). The number of available tiles in the following T est time slots can be estimated 
as T e (t) = J2n=i Et=o «™( t + r )> wnere ®n(t + 0) = a n (t) and t min = min{T eS f - l,T GoP - 
(t mod T GoP )}. T e (t) may not be an integer, but it does not affect the outcome of GRD2. 

We then adjust the Z ff , m 's based on T e {t) and N ac k(t), the number of ACKs received in time slot 
t. If T e {t) + N ac k(t — 1) > T e (t— 1)+N ac k(t— 2), there are more tiles that can be allocated and we 
can increase some of the / 9i?n 's. On the other hand, if T e (t) + N ack (t— 1) < T e (t— 1) + N ack (t — 2), 
we have to reduce some of the Z g , m 's. Due to layered videos, when we increase the number of 
allocated tiles, we only need to consider l g >m for m — rn', m! + 1, • • • , M, where MC m i is the 
highest MC scheme used in the previous time slot. Similarly, when we reduce the number of 
allocated tiles, we only need to consider l gim for m = m', rn! + 1, • • • , M. 

The refined greedy algorithm is given in Table Unl For time slot t, the complexity of GRD2 
is O(MGK), where K = \N ack (t - 1) - N ack (t - 2) + T e (t) - T e (t - 1)|. Since K < T e , the 
complexity of GRD2 is much lower than GRD1, suitable for execution in each time slot. 

3) Tile Scheduling in a Time Slot : In each time slot t, we need to schedule the remaining 
tiles for transmission on the N channels. We define Inc(g, m, i) to be the increase in the group 
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TABLE III 

The Refined Greedy Algorithm (GRD2) for Each Time Slot 



1: Initialize Z g , m = for all g and m 

2: Initialize A = {1,2, ••• ,G] 

3: Initialize N ack (0) = 

4: Estimate T e (l) based on the Markov Chain channel model 

5: Use GRD1 to find all l g , m 's based on T e (l) 

6: WHILE t = 2 to T GoP 

7: Estimate T e (t) 

8: IF [T e (i) + N ack (t - 1) < T e (t - 1) + W acfc (t - 2)] 

9: WHILE [£° =1 E™=i > T e {t) + ^a«*(* - 2)] 

10: Find that can be reduced by 1: 

_ • f U(T)-U(T-e„ m) 1 

e§,m - argmin VSim€{m ,,., M) \ <> S ,™+-R/Te j 

11. / — / £g,rh 

12: IF (3 A) 

13: Add g to A 

14: END IF 

15: END WHILE 

16: END IF 

17: IF [T e (t) + N ack (t - 1) > T e (i - 1) + N ack (t - 2)] 

18: WHILE [E° =1 E™=i < T e (t) + iV acfc (* - 1) and 

A is not empty] 

19: Find that can be increased by 1 

_ / U(l+S g m )-U(l) \ 
e g,rh — &rg Hax s gx,me{m',"-,M} 1 i> 3 m +-R/T e J 

20: f = f+ e§, A 

21: IF(E m ^,^, m >^) 

22. / — / ^g,m 

23: Delete <? from A 

24: END IF 

25: END WHILE 

26: END IF 

27: Update N ack (t - 1) 

28: END WHILE 



utility function U (g) after the i-th tile in the sub-layer using MC m is successfully decoded. It 
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TABLE IV 

Algorithm for Tile Scheduling in a Time Slot 


1: 


Initialize m g to the lowest MC that has not been ACKed for all g 


2: 


Initialize i g to the first packet that has not been ACKed for all g 


3: 


Sort {c„(t)} in decreasing order. Let the sorted channel list be 




indexed by j. 


4: 


While (j = 1 to N) 


5: 


Find the group having the maximum increase in U(g): 




g = arg maxv s lnc(g,rn g ,i g ) 


6: 


Allocate the tile on channel j to group g 


7: 


Update m g and i g 


8: 


End while 



can be shown that 

figbg,m 

Q^g ~t~ fig X^u=l ^9,u^g,u (j 1 tyPgbg,m 

Lnc(<7, m, i) can be interpreted as the reward if the tile is successfully received. 

Letting c n (t) be the probability that the tile is successfully received, then we have c n (t) = 
Pn(t)a n (t). Our objective of tile scheduling is to maximize the expected reward, i.e., 

N 

maximize: E[Reward(£)] = c n (t) ■ Inc(£ n ), (16) 

71=1 

where £ = {£ n }n=i,—,JV an d i s tne tile allocation for channel n, i.e., representing the three- 
tuple {g,m,i}. The TSA algorithm is shown in Table HVl which solves the above optimization 
problem. The complexity of TSA is O(NlogN). We have the following theorem for TSA. 
Theorem 2: E[Reward] is maximized if Inc(£j) > Inc(£j) when c$(t) > Cj(t) for all i and j. 
Proof: Suppose there exists a pair of i and j where Inc(£j) > Inc(£j) and Ci(t) < Cj(t). We 
can further increase E[Reward] by switching the tile assignment, i.e., assign channel i to £j and 
channel j to £j. With this new assignment, the net increase in E[Reward] is 

c i (t)Inc(£ i ) + Cj(i)Inc(£j) - c i (t)Inc(£ i ) - c i (t)Inc(£ i ) 
= [ Cj (t) - Ci {t)][lnc{Ci) - Inc(^)] > 0. 
Therefore E[Reward] is maximized when the {Inc(£j)} and {cj(t)} are in the same order. ■ 



M 

lnc(g,m,i) = ^{n gM - n fljfc+ i) x log 

k=m 
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C. Simulation Results 

We evaluate the performance of the proposed CR video multicast framework using a cus- 
tomized simulator implemented with a combination of C and MATLAB. Specifically, the LPs 
are solved using the MATLAB Optimization Toolbox and the remaining parts are written in C. 
For the results reported in this section, we have N = 12 channels (unless otherwise specified). 
The channel parameters A n and fx n are set between (0,1). The maximum allowed collision 
probability 7 n is set to 0.2 for all the channels unless otherwise specified. 

The CR base station multicasts three Common Intermediate Format (CIF, 352 x 288) video 
sequences to three multicast groups, i.e., Bus to group 1, Foreman to group 2, and Mother & 
Daughter to group 3. The n 1?n 's are {42, 40, 36, 30, 22, 12} (i.e., 42 users can decode MC% 
signal, 40 users can decode MCi signal, and so forth); the n 2) m , s are {51, 46, 40, 32, 23, 12} 
and the n 3 m 's are {49, 44, 40, 32, 24, 13}. The number of bits carried in one tile using the 
MC schemes are 1 kb/s, 1.5 kb/s, 2 kb/s, 3 kb/s, 5.3 kb/s, and 6 kb/s, respectively. We choose 
2c?oP= 150 and T est = 10, sensing interval W — 3, false alarm probability e n = 0.3 and miss 
detection probability S n = 0.25 for all n, unless otherwise specified. 

In every simulation, we compare three schemes: (i) a simple heuristic scheme that equally 
allocates tiles to each group (Equal Allocation); (ii) A scheme based on SF (Sequential Fixing), 
and (iii) a scheme based on the greedy algorithm GRD2 (Greedy Algorithm). These schemes 
have increasing complexity in the order of Equal Allocation, Greedy Algorithm, and Sequential 
Fixing. They differ on how to solve Problem OPT- Part, while the same tile scheduling algorithm 
and opportunistic spectrum access scheme are used in all the schemes. Each point in the figures 
is the average of 10 simulation runs, with 95% confidence intervals plotted. We observe that 
the 95% confidence intervals for Equal Allocation and Greedy Algorithm are negligible, while 
the 95% confidence intervals for Sequential Fixing is relatively larger. The C/MATLAB code is 
executed in a Dell Precision Workstation 390 with an Intel Core 2 Duo E6300 CPU working at 
1.86 GHz and a 1066 MB memory. For number of channels ranging from iV=3 to iV=15, the 
execution times of Equal Allocation and Greedy Algorithm are about a few milliseconds, while 
Sequential Fixing takes about two seconds. 

In Fig. [5] we plot the average PSNR among all users in each multicast group. For all the 
groups, Greedy Algorithm achieves the best performance, with up to 4.2 dB improvements over 
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Group Index 



Fig. 5. Average PSNR of all multicast users. 

Equal Allocation and up to 0.6 dB improvements over Sequential Fixing. We find Sequential 
Fixing achieves a lower PSNR than Equal Allocation for group 3, but higher PSNRs for groups 
1 and 2. This is because Equal Allocation does not consider channel conditions and fairness. 
It achieves better performance for group 3 at the cost of much lower PSNRs for groups 1 and 
2. We also plot Frame 53 from the original Bus sequence and the decoded video at user 1 of 
group 1 in Fig [6l We choose this user since it is one of the users with lowest PSNR values. The 
average PSNR of this user is 29.54 dB, while the average PSNR of all group 1 users is 34.6 dB. 
Compared to the original frame (right), the reconstructed frame (left) looks quite good, although 
some details are lost. 

In Fig. |71 we examine the impact of the maximum allowed collision probability 7 n . We increase 
7 n from 0.1 to 0.3, and plot the average PSNR values among all the users. When 7„ gets larger, 
there will be higher chance of collision for the video packets, which hurts the received video 
quality. However, a higher 7 n also allows a higher transmission probability Pn(t) for the base 
station (see ©), thus allowing the base station to grab more spectrum opportunities and achieve 
a higher video rate. The net effect of these two contradicting effects is improved video quality 
for the range of 7 n values considered in this simulation. This is illustrated in the figure where all 
the three curves increase as 7 n gets larger. We also observe that the curves for Sequential Fixing 
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Fig. 6. The original (the right one) and decoded Frame 53 (the left one) at user 1 in group 1. 
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Fig. 7. Average PSNR of all users versus 7„ (with 95% confidence intervals). 

and Equal Allocation are roughly parallel to each other, while the Greedy Algorithm curve has a 
steeper slope. This indicates that Greedy Algorithm is more efficient in exploiting the additional 
bandwidth allowed by an increased j n . 

In Fig. [8} we examine the impact of number of channels N. We increase N from 3 to 15 
in steps of 3, and plot the average PSNR values of all multicast users. As expected, the more 
channels, the more spectrum opportunities for the CR networks, and the better the video quality. 
Again, we observe that the Greedy Algorithm curve has the steepest slope, implying it is more 
efficient in exploiting the increased spectrum opportunity for video transmissions. 
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Fig. 8. Average PSNR of all users versus N (with 95% confidence intervals). 



We demonstrate the impact of sensing errors in Fig. [9j We test five sets of {e n ,S n } values as 
follows: {0.10,0.38}, {0.30,0.25}, {0.5,0.17}, {0.70,0.10} and {0.9,0.04} 02, and plot the 
average PSNR values of all users. It is quite interesting to see that the video quality is not very 
sensitive to sensing errors. Even as t n is increased nine times from 10% to 90%, there is only 
0.58 dB reduction (or a 1.5% normalized reduction) in average PSNR when Greedy Algorithm 
is used. The same can be observed for the other two curves. We conjecture that this is due to the 
opportunistic spectrum access approach adopted in all the three schemes. A special strength of 
the proposed approach is that it explicitly considers both types of sensing errors and mitigates 
the impact of both sensing errors. For example, when the false alarm rate is very high, the base 
station will not trust the sensing results and will access the channel relatively more aggressively, 
thus mitigating the negative effect of the high false alarm rate. 

Finally, we demonstrate the impact of user channel variations (i.e., due to mobility). We chose 
a tagged user in group 1 and assume that its channel condition changes every 20 GoPs. The 
highest MC scheme that the tagged user can decode is changed according to the following 
sequence: MC3, MC5, MC4, MC6, MC5 and MC3. All other parameters remain the same as in 
the previous experiments. In Fig. [TQl we plot the average PSNRs for each GoP at this user that 
are obtained using the three algorithms. We observe that both Greedy Algorithm and Sequential 
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Fig. 9. Average PSNR of all users for various {e n ,5 n } values (with 95% confidence intervals). 



Fixing can quickly adapt to changing channel conditions. Both algorithms achieve received video 
qualities commensurate with the channel quality of the tagged user. We also find the video quality 
achieved by Greedy Algorithm is more stable than that of Sequential Fixing, while the latter 
curve has some deep fades from time to time. This is due to the fact that Greedy Algorithm 
has a proven optimality bound, while Sequential Fixing does not provide any guarantee. The 
Equal Allocation curve is relative constant for the entire period since it does not adapt to channel 
variations. Although being simple, it does not provide good video quality in this case. 

For optimization-driven multimedia systems, there is a trade-off between (i) grabbing all the 
available resource to maximize media quality and (ii) be less adaptive to network dynamics for 
a smooth playout. The main objective of this paper is to demonstrate the feasibility and layout 
the framework for video streaming over infrastructure-based CR networks, using an objective 
function of maximizing the overall user utility. We will investigate the interesting problem of 
trading off resource utilization and smoothness in our future work. 

V. Video over Multi-hop CR Networks 

In this section, we examine the problem of video over multi-hop CR networks. We model 
streaming of concurrent videos as an MINLP problem, aiming to maximize the overall received 



24 



44 




42- 



Greedy Algorithm 

Sequential Fixing 

Equal Allocation 



28 







20 



40 



60 

GoP Index 



80 



100 



120 



Fig. 10. GoP average PSNRs of a tagged user in Group 1, when its channel condition varies over time. 

video quality and fairness among the video sessions, while bound the collision rate with primary 
users under spectrum sensing errors. We solve the MINLP problem using a centralized sequential 
fixing algorithm, and derive upper and lower bounds for the objective value. We then apply 
dual decomposition to develop a distributed algorithm and prove its optimality and convergence 
conditions. 

A. Network Model 

1) Spectrum Access : During the transmission phase of a time slot, a CR user determines 
which channel(s) to access for transmission of video data based on spectrum sensing results. 
Let K k m be a threshold for spectrum access: channel m is considered idle if the estimate is 
greater than the threshold, and busy otherwise. The availability of channel m in primary network 
k, denoted as At,, is 



For each channel m, we can calculate the probability of collision with primary users as: 




1, otherwise. 



(17) 




(18) 
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where set ^ is defined as: 

%b k = 



> k1 



(19) 



Pr(^o|< 

For non-intrusive spectrum access, the collision probability should be bounded with a prescribed 
threshold 7^. A higher spectrum access threshold n k m will reduce the potential interference with 
primary users, but increase the chance of wasting transmission opportunities. For a given collision 
tolerance 7^, we can solve Pr(A^ = 0\Hi) = 7^ for n k m . The objective is to maximize CR 
users' spectrum access without exceeding the maximum collision probability with primary users. 
Let be the set of available channels at link {i,j}. Assuming i eU k and j G U k ' , we have 



a 



{m |^ = and ^ = 0}. 



(20) 



2) Link and Path Statistics: Due to the amplify-and-forward approach for video data trans- 
mission, there is no queueing delay at intermediate nodes. Assume each link has a fixed delay 
Uij (i.e., processing and propagation delays). Let Vf be the set of all possible paths from z\ 
to d\. For a given delay requirement T th , the set of feasible paths V\ for video session I can be 
determined as: 



V 



E 



< T th , V e V{ 



(21) 



Let p™j be the packet loss rate on channel m at link {i,j}. A packet is successfully delivered 
over link {i,j} if there is no loss on all the channels that were used for transmitting the packet. 
The link loss probability p i:j can be derived as: 



P 



=i - n (i-co^ 



(22) 



where M. is set of licensed channels and I rn is an indicator: I m — 1 if channel m is used for 
the transmission, and I m — otherwise. Assuming independent link losses, the end-to-end loss 
probability for path Vj 1 G V\ can be estimated as: 



(23) 
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B. Problem Statement 

We also aim to achieve fairness among the concurrent video sessions. It has been shown that 
proportional fairness can be achieved by maximizing the sum of logarithms of video PSNRs 
(i.e., utilities). Therefore, our objective is to maximize the overall system utility, i.e., 

maximize: ^ U t (Ri) = ^log(Q,0Ri))- (24) 
i i 

1) Multi-hop CR Network Video Streaming Problem: The problem of video over multi-hop 
CR networks consists of path selection for each video session and channel scheduling for each 
CR node along the chosen paths. We define two sets of index variables. For channel scheduling, 
we have 

1, at link {i, j}, if channel m is 

assigned to tunnel r in path Vj 1 (25) 
0, otherwise. 



Lh.r 



For path selection, we have 

, 1, if video session / selects path Vj 1 G Vi 

!/?=< ' (26) 

! 0, otherwise, 

Note that the indicators, x^ T m and y\, are not independent. If yf = for path V[\ all the 
xjj'^'s on that path are 0. If link {i, j} is not on path Vj 1 , all its a^'^'s are also 0. For link 
{i,j} on path P/ 1 , we can only choose those available channels in set fijj to schedule video 
transmission. That is, we have x^'^ G {0, 1} if m G and x 1 -^ = otherwise. In the rest 
of the paper, we use x and y to represent the vector forms of x^'^ and respectively. 

As discussed, the objective is to maximize the expected utility sum at the end of A^ G time 
slots, as given in (|24l) . Since \og(Qi(E[Ri(0)])) is a constant, d24l) is equivalent to the sum of 
utility increments of all the time slots, as 



5>g(Q,(E[J?,(JV )])) - tog(g,(E[i?,(0)])) 
i 

E {log(QKE[^(*)D) - log(Q,(E[J2,(i - 1)]))} • (27) 
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Therefore, (1241) will be maximized if we maximize the expected utility increment during each 
time slot, which can be written as: 



E log(Q,(E[J2,(f)])) - log(0*(B[J%(* - 1)])) 

^log(i + A E ^ (t)] - E[ ^ (t - 1)r 



Qi(E[Ri(t-l)]) 



l.h.r 



NgTsQi 



E E f? lo s f 1+ ^ E E ^.J 1 - pw) 



where Zj is the next hop from z\ on path V\\ p\ h is the packet loss rate on tunnel r of path Vj 1 , 
Q\- 1 = Qi^R^t - 1)]), and p\ = ^/{NgTsQ^ 1 ). 

From ((22)) and (1231) . the end-to-end packet loss rate for tunnel r on path Vj 1 is: 

Pi,ft = 1 - II H^-PZT*^- (28) 

{i,i}G7 , i ?l mGA4 

We assume that each tunnel can only include one channel on each link. When there are multiple 
channels available at each link along the path, a CR source node can set up multiple tunnels to 
exploit the additional bandwidth. We then have the following constraint: 

X>&»<1. V{z,j}e7t (29) 

Considering availability of the channels, we further have, 

EE41m<l^,il, V{iJ}eV?, (30) 

r m 

where is the number of available channels on link {i,j} defined in (|2Q|) . 

As discussed, each node is equipped with two transceivers: one for receiving and the other 
for transmitting video data during the transmission phase. Hence a channel cannot be used to 
receive and transmit data simultaneously at a relay node. We have for each channel m: 

E *ftm + E ■'■'£»■ ^ X > V m, /, V G P,, V {z, j}, {j, k} G Vl (31) 



28 



Let nf be the number of tunnels on path Vj 1 . For each source z\ and each destination di, the 
number of scheduled channels is equal to nf. We have for each source node 

EE<X™ = ^' VhEP h Vl. (32) 

r m 

Let d\ be the last hop to destination di on path Vj 1 , we have for each destination node 

EE^J,^^ Phenyl. (33) 

r m 

At a relay node, the number of channels used to receive data is equal to that of channels used 
to transmit data, due to flow conservation and amplify-and-forward. At relay node j for session 
I, assume {i, j} G Vj 1 and {j, k} G Vf. We have, 

E E ^ = E E 4ii> v /i g v /, v o, j}, u fc} g t?. (34) 

r m r m 

We also consider hardware-related constraints on path selection. We summarize such con- 
straints in the following general form for ease of presentation: 

To simplify exposition, we choose at most one path in Vi for video session I. Such a single 
path routing constraint can be expressed as YlhVl — 1> wn i cn i s a special case of (|35l) where 
= 1 for all h, and wf, fe = for all g / 1, I' 7^ i, and h. We can also have YlihVi ^ ^ t0 
allow up to £ paths for each video session. In order to achieve optimality in the general case of 
multi-path routing, an optimal scheduling algorithm should be designed to dispatch packets to 
paths with different conditions (e.g., different number of tunnels and delays). 

There are also disjointedness constraints for the chosen paths. This is because each CR node 
is equipped with two transceivers and both will be used for a video session if it is included in 
a chosen path. Such disjointedness constraint is also a special case of (T331) with the following 
definition for wf h for each CR node g: 

1, if node g G path V\ l 
0, otherwise 

Finally we formulate the problem of multi-hop CR network video streaming (OPT-CRV) as: 
max: y\ log (l + p\ £ £ x l ^Jl -p[ h )\ (37) 

l h&Vi \ r m J 

subject to: ([25]) ~ (135]) . 



< h ={ 1 (36) 
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2) Centralized Algorithm and Upper/Lower Bounds: Problem OPT-CRV is in the form of 
MINLP (without continuous variables), which is NP-hard in general. We first describe a cen- 
tralized algorithm to derive performance bounds in this section, and then present a distributed 
algorithm based on dual decomposition in the next section. 

We first obtain a relaxed non-linear programming (NLP) version of OPT-CRV. The binary 
variables x^/ m and y\ are relaxed to take values in [0,1]. The integer variables nf are treated 
as nonnegative real numbers. It can be shown that the relaxed problem has a concave object 
function and the constraints are convex. This relaxed problem can be solved using a constrained 
nonlinear optimization problem solver. If all the variables are integer in the solution, then we 
have the exact optimal solution. Otherwise, we obtain an infeasible solution, which produces an 
upper bound for the problem. This is given in Lines 1 ~2 in Table |V] 

We also develop a sequential fixing algorithm (SF) for solving OPT-CRV. The pseudo-code is 
given in Table |V] SF iteratively solves the relaxed problem, fixing one or more integer variables 
after each iteration [26]. In Table fVl Lines 3~7 fix the path selection variables yj 1 , and Lines 
8~16 fix the channel scheduling variables x 1 -^ and tunnel variables nf. The tunnel variables 
n\ can be computed using (1321) after x l ^ m and y\ are solved. When the algorithm terminates, 
it produces a feasible solution that yields a lower bound for the objective value. 

C. Dual Decomposition 

SF is a centralized algorithm requiring global information. It may not be suitable for multi- 
hop wireless networks, although the upper and lower bounds provide useful insights on the 
performance limits. In this section, we develop a distributed algorithm for Problem OPT-CRV 
and analyze its optimality and convergence performance. 

1) Decompose Problem OPT-CRV: Since the domains of x\ h ^ m defined in (|29l) ~(l34l) for 
different paths do not intersect with each other, we can decompose Problem OPT-CRV into two 
subproblems. The first subproblem deals with channel scheduling for maximizing the expected 
utility on a chosen path V\ l . We have the channel scheduling problem (OPT-CS) as: 




(38) 



r 



in 



subject to: plj ~ (JSJ, x 



l,h,r 
Z[,z' l; m 



G {0, 1}, for all /, h, r, m. 
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TABLE V 

The Sequential Fixing Algorithm (SF) for Problem OPT-CRV 



1 : 


Relax integer variables x^'^, yf, and nf; 


2 : 


Solve the relaxed problem using a constrained NLP solver; 


3 : 


if (there is y\ not fixed) 


4 : 


Find the largest y v , where [l',fi] = argmaxji/f}, 




and fix it to 1; 


5 : 


Fix other yf's according to constraint d35t: 


6 : 


Go to Step 2; 


7 : 


end if 


8 : 


if (there is x\^ m not fixed) 


9 : 


Find the largest x\,' ,' T m i, where m', I', h', r'] = 




argmax{x''^}, and set it to 1; 


10: 


Fix other x^'^'s according to the constraints; 


11: 


if (there is other variable that is not fixed) 


12: 


Go to Step 2; 


13: 


else 


14: 


Fix nj"s based on x and y; 


15: 


Exit with feasible solution {x, y,n}; 


16: 


end if 


17: 


end if 



In the second part, optimal paths are selected to maximize the overall objective function. Letting 
F[ l = log (l + pfHj 1 ), we have the following path selection problem (OPT-PS): 

maximize: /(y) = J2J2 < 39 ) 

l h 

subject to: w ihVi' — ^' f° r a ^ 9 

i heVi 

y\ E {0,1}, for all l,h. 

2) Solve the Channel Scheduling Subproblem: We have the following result for assigning 
available channels at a relay node. 

Theorem 3: Consider three consecutive nodes along a path, denoted as nodes i, j, and k. Idle 
channels 1 and 2 are available at link {i, j} and idle channels 3 and 4 are available at link {j, k}. 
Assume the packet loss rates of the four channels satisfy pjj > p?j and p^ k > pj k . To set up 
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two tunnels, assigning channels {1, 3} to one tunnel and channels {2, 4} to the other tunnel 
achieves the maximum expectation of successful transmission on path section k}. 

Proof: Let the success probabilities on the channels be pjj = 1 — p\ j, p\j = 1 — p? •, 

P\k — 1 ~~ P 3 j,k' an d Pj t k = 1 _ Pj,k- We nave Pi,j < P%j an ^ P%h < Pj,k- Comparing the success 
probabilities of the channel assignment given in Theorem |3] and that of the alternative assignment, 

we have Pi d 0j, k +PijPj,k-PijPj,k ~PijPj,k = {Pi,j-Pi,j)(plk -Pj,k) > °- The result follows. ■ 
According to Theorem [3j a greedy approach, which always chooses the channel with the 
lowest loss rate at each link when setting up tunnels along a path, produces the optimal overall 
success probability. More specifically, when there is only one tunnel to be set up along a path, 
the tunnel should consist of the most reliable channels available at each link along the path. 
When there are multiple tunnels to set up along a path, tunnel 1 should consist of the most 
reliable channels that are available at each link; tunnel 2 should consist of the second most 
reliable links available at each link; and so forth. 

Define the set of loss rates of the available channels on link {i, j} as Aj j = {p™j\m G f^}. 
The greedy algorithm is given in Table |VH with which each video source node solves Problem 
OPT-CS for each feasible path. Lines 2~3 in Table IVTl checks if there is more channels to assign 
and the algorithm terminates if no channel is left. In Lines 4~10, links with only one available 
channel are assigned to tunnel r and the neighboring links with the same available channels are 
removed due to constraint (f3TT) . In Lines 1 1 ~ 17, links with more than two channels are grouped 
to be assigned later. In Lines 18~20, the available channel with the lowest packet loss rate is 
assigned to tunnel r at each unallocated link, according to Theorem |3l To avoid co-channel 
interference, the same channel on neighboring links is removed as in Lines 21^33. 

3) Solve the Path Selection Subproblem: To solve Problem OPT- PS, we first relax binary 
variables y\ to allow them take real values in [0,1] and obtain the following relaxed path 
selection problem (OPT-rPS): 

maximize: /(y) = J2Y1 F iVi < 40) 

/ h 

subject to: w fhVi — 1' f° r a ^ 9 

i heVi 

0<y?< 1, for all h, I 
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We then introduce positive Lagrange Multipliers e g for the path selection constraints in Problem 
OPT-rPS and obtain the corresponding Lagrangian function: 

£<y> e ) = E E F M + E e ^ - E E <yi) ( 41 ) 

l h g I h 

= EE^v-E<^) + E e ^ 

l h g g 

l h g 

Problem (|4TT) can be decoupled since the domains of y^s do not overlap. Relaxing the coupling 
constraints, it can be decomposed into two levels. At the lower level, we have the following 
subproblems, one for each path Vj 1 , 



max C^vt e) = f*y* - £ jfo. (42) 
At the higher level, by updating the dual variables e g , we can solve the relaxed dual problem: 

mm g(e) = E E £ " ( MT > e ) + E < 43 > 

i ft g 

where is the optimal solution to (|42)) . Since the solution to (142)) is unique, the relaxed dual 

problem (|43T) can be solved using the following subgradient method that iteratively updates the 
Lagrange Multipliers [|29ll : 

1 + 

^~ ' , (44) 



e (r + L 



e g (r)-a(r)(l-J2J2 w Uyi) 

l h 

where r is the iteration index, a(r) is a sufficiently small positive step size and [x] + denotes 
max{x, 0}. The pseudo code for the distributed algorithm is given in Table rVTfl 

4) Optimality and Convergence Analysis: The distributed algorithm in Table [VTll iteratively 
updates the dual variables until they converge to stable values. In this section, we first prove that 
the solution obtained by the distributed algorithm is also optimal for the original path selection 
problem OPT-PS. We then derive the convergence condition for the distributed algorithm. 

Fact 1 ( K29\l): Consider a linear problem involving both equality and inequality constraints 

maximize: a'x (45) 
subject to: h^x = b%, ■ ■ ■ , h^x = b m 
g'lX < d, • • •, g',x < c r , 
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where a, hj, and gj are column vectors in lZ n , frj's and c/s are scalars, and a' is the transpose 
of a. For any feasible point x, the set of active inequality constraints is denoted by *4(x) = 
{j|g^x = Cj}. If x* is a maximizer of inequality constrained problem (1451 ), x* is also a maximizer 
of the following equality constrained problem: 

maximize: a'x (46) 

subject to: h^x = b\, ■ ■ ■ , h^x = b rn 

g;x = Cj ,vj eA(x). 

Lemma 1: The optimal solution for the relaxed primal problem OPT-rPS in (|4Q|) is also feasible 
and optimal for the original Problem OPT- PS in (|39l) . 

Proof: According to Fact [fl the linearized problem of OPT-PS, i.e., OPT-rPS, can be 
rewritten as an equality constrained problem in the following form: 

maximize: F'y (47) 

subject to: w^-y = 1, j G A(y*) (48) 

< Vi < 1, for all h,l, 

where F, w/s, and y are column vectors with elements F t h , wf h , and y\, respectively. We apply 
Gauss-Jordan elimination to the constraints in (l48l to solve for y. Since there is not sufficient 
number of equations, some y\ 's are free variables (denoted as y{) and the rest are dependent 
variables (denoted as yf). Assuming there are r free variables, the dependent variables can be 
written as linear combinations of the free variables after Gauss-Jordan elimination, as 

r 
i=l 

Due to Gauss-Jordan elimination and binary vectors w/s, Wj and bj in (|49l) are all integers. 
Therefore, if all the free variables y{ attain binary values, then all the dependent variables 
computed using (|49| ) will also be integers. Since < yj < 1, being integers means that they are 
either or 1, i.e., binaries. That is, such a solution will be feasible. 

Next we substitute (l49l) into problem (l47l) to eliminate all the dependent variables. Then we 
obtain a unconstrained problem with only r free variables, as 

r 

maximize: 2J Fiy( + b (50) 

i=l 
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Since the free variables yf's take value in {0, 1}, this problem can be easily solved as follows. If 
the coefficient Fj > 0, we set y{ = 1; otherwise, if Fj < 0, we set y{ = 0. Thus (l50l) achieves its 
maximum objective value. Once all the free variables are determined with their optimal binary 
values, we computes the dependent variables using (|49|) . which are also binary as discussed 
above. Thus we obtain a feasible solution, which is optimal. ■ 

Lemma 2: If the relaxed primal Problem OPT-rPS in (|40l) has an optimal solution, then the 
relaxed dual problem (1431) also has an optimal solution and the corresponding optimal values of 
the two problems are identical. 

Proof: By definition, the problems in (14TT) and (1431) are primal/dual problems. The primal 
problem always has an optimal solution because it is bounded. Since Problem OPT-rPS is an 
LP problem, the relaxed dual problem is also bounded and feasible. Therefore the relaxed dual 
problem also has an optimal solution. We have the strong duality if the primal problem is convex, 
which is the case here since Problem OPT-rPS is an LP problem. ■ 

We have Theorem 0] on the optimality of the path selection solution, which follows naturally 
from Lemmas Q] and [2l 

Theorem 4: The optimal solution to the relaxed dual problem (|42l) and (|43l) is also feasible 
and optimal to the original path selection Problem OPT- PS given in (|39| ). 

As discussed, the relaxed dual problem (|43l) can be solved using the subgradient method that 
iteratively updates the Lagrange Multipliers. We have the following theorem on the convergence 
of the distributed algorithm given in Table IVIII 

Theorem 5: Let e* be the optimal solution. The distributed algorithm in Table IVIII converges 
if the step sizes a(r) in (l44l) satisfy the following condition: 

< a(r) < ^'tt'Tnf" " ■ for all r, (51) 

\\G{t)\\ 2 

where G(r) is the gradient of g(e(r)). 

Proof: Since g(e(r)) is a linear function, we have subgradient equality, as 

g (e(r))-g(e*) = [e(r)-e*] / G(r). 
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It then follows that 

||e(r) -a(r)G(r) -e*|| 2 

= ||e(r)-e*|| 2 - 2a(r)[e(r) -e*]'G(r) + (a(r)) 2 ||G(r)|| 2 
= ||e(r)-e*|| 2 - 2a(r)[g(e(r))-g(e*)] + (a(r)) 2 ||G(r)|| 2 

(52) 

If a(r) satisfy (ISTI) . the sum of the last two terms in (|52|) is negative. It follows that, ||e(r) — 
a(r)G(r) — e*|| < ||e(r) — e*||. Since the projection operation is nonexpansive, we have, 

| |e(r + 1) - e* 1 1 = 1 1 [e(r) - a(r)G(r)] + - [e*]+| | 
< |)e(r) - a(r)G(r) - e*|| < ||e(r) - e*||, 

which states the conditional convergence of the algorithm. ■ 
Since the optimal solution e* is not known a priori, we use the following approximation in 
the algorithm: a{r) = ^fclp, where q(r) is the current estimate for q(e*). We choose the 
mean of the objective values of the relaxed primal and dual problems for q(r). 

5) Practical Considerations: Our distributed algorithms are based on the fact that the com- 
putation is distributed on each feasible path. The OPT-CS algorithm requires information on 
channel availability and packet loss rates at the links of feasible paths. The OPT-PS algorithm 
computes the primal variable y\ for each path and broadcasts Lagrangian multipliers over the 
control channel to all the source nodes. We assume a perfect control channel such that channel 
information can be effectively distributed and shared, which is not confined by the time slot 
structure 03). 

We assume relatively large timescales for the primary network time slots, and small to medium 
diameter for the CR network, such that there is sufficient time for timely feedback of channel 
information to the video source nodes and for the distributed algorithms to converge. Otherwise, 
channel information can be estimated using © based on delayed feedback, leading to suboptimal 
solutions. If the time slot is too short, the distributed algorithm may not converge to the optimal 
solution (see Fig. [15]). We focus on developing the CR video framework in this paper, and will 
investigate these issues in our future work. 
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Fig. 11. Topology of the multi-hop CR network. Note that only video source nodes, video destination nodes, and those nodes 
along the precomputed paths are shown in the topology. 



D. Simulation Results 

1) Methodology and Simulation Settings: We implement the proposed algorithms with a 
combination of C and MATLAB (i.e., for solving the relaxed NLP problems), and evaluate 
their performance with simulations. For the results reported in this section, we have K = 3 
primary networks and M — 10 channels. There are 56, 55, and 62 CR users in the coverage 
areas of primary networks 1, 2, and 3, respectively. The |W^J's are [5 46487567 4] (i.e., 
five users sense channel 1, four users sense channel 2, and so forth); the |W^J's are [4 6 5 7 6 
5 3 8 5 6], and the |Z^|'s are [8 6 5 4 7 6 8 5 6 7]. The topology is shown in Fig. Q3J 

We choose L p = 100, T s = 0.02 and N G = 10. The channel utilization is r\ k m = 0.6 for all 
the channels. The probability of false alarm is e^ n = 0.3 and the probability of miss detection 
is 5^ = 0.2 for all m and k, unless otherwise specified. Channel parameters and [i k m are 
set between (0, 1). The maximum allowed collision probability 7^ is set to 0.2 for all the M 
channels in the three primary networks. 

We consider three video sessions, each streaming a video in the Common Intermediate Format 
(CTF, 352 x 288), i.e., Bus to destination 1, Foreman to destination 2, and Mother & Daughter 
to destination 3. The frame rate is 30 fps, and a GOP consists of 10 frames. We assume that the 
duration of a time slot is 0.02 seconds and each GOP should be delivered in 0.2 seconds (i.e., 
10 time slots). 
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We compare four schemes in the simulations: (i) the upper-bounding solution by solving 
the relaxed version of Problem OPT-CRV using an NLP solver, (ii) the proposed distributed 
algorithm in Tables [VI] and IVIIl (iii) the sequential fixing algorithm given in Table [V] which 
computes a lower-bounding solution, and (iv) a greedy heuristic where at each hop, the link with 
the most available channels is used. Each point in the figures is the average of 10 simulation 
runs, with 95% confidence intervals plotted as error bars in the figures. The 95% confidence 
intervals are negligible in all the figures. 

2) Simulation Results: 

a) Algorithm Performance: To demonstrate the convergence of the distributed algorithm, we 
plot the traces of the four Lagrangian multipliers in Fig. [[2] We observe that all the Lagrangian 
multipliers converge to their optimal values after 76 iterations. We also plot the control overhead 
as measured by the number of distinct broadcast messages for e^r) using the y-axis on the right- 
hand side. The overhead curve increases linearly with the number of iterations and gets flat (i.e., 
no more broadcast message) when all the Lagrangian multipliers converge to their optimal values. 

We examine the impact of spectrum sensing errors in Fig. [13] We test six sensing error 
combinations {e m ,5 m } as follows: {0.1, 0.5}, {0.2, 0.3}, {0.3, 0.2}, {0.5, 0.11}, {0.7, 0.06}, 
and {0.9, 0.02}, and plot the average PSNR values of the Foreman session. It is interesting to 
see that the best video quality is achieved when the false alarm probability e m is between 0.2 
and 0.3. Since the two error probabilities are correlated, increasing one will generally decrease 
the other. With a larger e m , CR users are more likely to waste spectrum opportunities that are 
actually available, leading to lower bandwidth for videos and poorer video quality, as shown 
in Fig. [13] On the other hand, a larger 5 m implies more aggressive spectrum access and more 
severe interference to primary users. Therefore when e m is lower than 0.2 (and 5 rn is higher than 
0.3), the CR nodes themselves also suffer from the collisions and the video quality degrades. 

b) Impact of Primary Network Parameters: In Fig. [Q] we examine the impact of channel 
utilization rj on received video quality. We focus on Session 2 with the Foreman sequence. The 
average PSNRs achieved by the four schemes are plotted when 77 is increased from 0.6 to 0.9 
for all licensed channels. Intuitively, a smaller i] allows more transmission opportunities for CR 
nodes, leading to improved video quality. This is illustrated in the figure where all the four curves 
decrease as 77 gets larger. The distributed scheme achieves PSNRs very close to that obtained 
by sequential fixing, and both of them are close to the upper bound. The heuristic scheme is 
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Fig. 12. Illustrate the convergence of the distributed algorithm. 
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Fig. 13. Video PSNRs versus spectrum sensing error. 



inefficient in exploiting the available spectrum even when the channel utilization is low. As 
discussed, the time slot duration is also an important parameter that may affect the convergence 
of the distributed algorithm. In Fig. \\5i we keep the same network and video session settings, 
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Fig. 14. Video PSNRs versus primary user channel utilization r\. 



while increasing the time slot duration as 4 ms, 10 ms, 20 ms, 40ms and 100 ms. For a given 
time slot duration, we let the distributed algorithm run for 5% of the time slot duration, starting 
from the beginning of the time slot, and then stop. The solution that the algorithm produces 
when it is stopped will be used for video transmission in the remainder of this time slot. It can 
be seen that when the time slot is 4 ms, the algorithm does not converge after 5% x 4=0.2 ms, and 
the PSNR produced by the distributed algorithm is low (but still higher than that of the heuristic 
algorithm). When the time slot duration is sufficiently large (e.g., over 10 ms), the algorithm 
can converge and the proposed algorithm produces very good video quality as compared to the 
upper bound and the lower bound given by the sequential fixing algorithm. 

c) Comparison ofMPEG-4 FGS and H.264/SVC MGS Videos: Finally, we compare MPEG- 
4 FGS and H.264/SVC MGS videos, while keeping the same settings. It has been shown that 
H.264/SVC has better rate-distortion performance than MPEG-4 FGS due to the use of efficient 
hierarchical prediction structures, the inter-layer prediction mechanisms, improved drift control 
mechanism, and the efficient coding scheme in H.264/AVC 0. Although MGS has Network 
Abstraction Layer (NAL) unit-based granularity, it achieves similar rate-distortion performance 
as H.264/SVC FGS 0. 

We plot the upper bounds and the distributed algorithm results in Figs. [J6] and [17] for vari- 
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Fig. 15. Impact of time slot duration on received video quality. 

ous channel utilizations and false alarm probabilities, respectively. From the figures, it can be 
observed that there is a gap about 2.5 dB between the H.264/SVC MGS and MPEG-4 FGS 
curves, which clearly demonstrates the rate-distortion efficiency of MGS over MPEG-4 FGS. 
The proposed algorithm can effectively handle both MGS and FGS videos, and the same trend 
is observed in both cases. 

VI. Conclusions 

In this paper, we first addressed the problem of multicasting FGS video in CR networks. The 
problem formulation took video quality and proportional fairness as objectives, while considering 
cross-layer design factors such as FGS coding, spectrum sensing, opportunistic spectrum access, 
primary user protection, scheduling, error control and modulation. We proposed efficient opti- 
mization and scheduling algorithms for highly competitive solutions, and proved the complexity 
and optimality bound of the proposed greedy algorithm. Our simulation results demonstrate not 
only the viability of video over CR networks, but also the efficacy of the proposed approach. 

Then, we studied the challenging problem of streaming multiple scalable videos in a multi-hop 
CR network. The problem formulation considered spectrum sensing and sensing errors, spectrum 
access and primary user protection, video quality and fairness, and channel/path selection for 
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Fig. 16. Comparison of MPEG-4 FGS video with H.264/SVC MGS video under various channel utilizations. 



48 
46 
44 
42 

ST 

S. 40 
CC 

CO 38 
0_ 

36 
34 
32 
30 



X— _- 


s - 


Upper Bound: MGS 
Distributed: MGS 


Upper Bound: FGS 


\ \ N 

/ \s N> 
\v 


Distributed: FGS 













0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
False Alarm Probability (e) 



Fig. 17. Comparison of MPEG-4 FGS video with H.264/SVC MGS video under various false alarm probabilities. 



concurrent video sessions. We first solved the formulated MINLP problem using a sequential 
fixing scheme that produces lower and upper bounds on the achievable video quality. We then 
applied dual decomposition to derive a distributed algorithm, and analyzed its optimality and 
convergence performance. Our simulations validated the efficacy of the proposed scheme. 
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TABLE VI 

The Greedy Algorithm for Channel Scheduling 



1 : Initialization: tunnel r = 1, link {i,j}'s from zi to di; 

2: if (lA^I ==0) 

3 : Exit; 

4: else if (\A t>j \ == 1) 

5 : Assign the single channel in A,,.,-, m' , to tunnel r; 

6 Check neighboring link {k, i}; 

7 : if (pj£; € A M ) 

8 : Remove from A fc>i , 

i <— fc, j <r- i and go to Step 2; 

9 : else 

10: Go to Step 13; 

11: end if 

12: else 

13: Put Aij in set Af; 

14: if (node j is not destination d t ) 

15: i <— j, j <— w; 

16: Go to Step 2; 

17: end if 

18: end if 

19: while (Af is not empty) 

20: Find the maximum value p™j, in set Af 

{i',j',m'} = argminjp™,}; 

21: Assign channel m' to tunnel r; 

22: Remove set A;/y from set Af; 

23: Check neighboring link {k, i} and {j, v}; 

24: if (pj$ e A fcli and A fc>< G Af) 

25: Remove pjl^j from Aft,;; 

26: if (Ak,i is empty) 

27: Exit; 

28: end if 

29: end if 

30: if ipf; v e Aj, v and A,,, € Af) 

31: Remove p™ v from Aj :V ; 

32: if (Aj ;V is empty) 

33: Exit; 

34: end if 

35: end if 

36: end while 

37: Compute the next tunnel: r «— r + 1 and go to Step 2; 
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TABLE VII 

Distribution Algorithm for Path Selection 



1: Initialization: set r = 0, e 3 (0) > and step size s G [0, 1]; 
2: Each source locally solves the lower level problem in J42t ; 

if (Ft - E 9 dl h e g {r)) > 0) yf = y? + a> y? = minfa* 1}; 
else yf = yf - s, y, h = max{yf , 0}; 
3: Broadcast solution y i ' l (e(T)); 

4: Each source updates e according to ( 1441 ) and broadcasts e(r + 1) 

through the common control channel; 
5: t <— t+1 and go to Step 2 until termination criterion is satisfied; 
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